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Abstract 

We show that any first order ordinary differential equation with a known Lie 
point symmetry group can be discretized into a difference scheme with the same 
symmetry group. In general, the lattices are not regular ones, but must be adapted 
to the symmetries considered. The invariant difference schemes can be so chosen 
that their solutions coincide exactly with those of the original differential equation. 



1 Introduction 

The purpose of this article is to analyze Lie point symmetries and to obtain exact solu- 
tions of first order difference schemes. These are two-point schemes of the form: 

Ea{x,y,x+,y+) ^ 0, a =1,2, (1.1) 

d(Ei,E2) , 
^ ' ' 7^0, (1.2) 



5(a;+,y+) 

where we use the notation x = Xn, x+ = x„+i, y = yn, y+ = 2/ri+i for the indepen- 
dent and dependent variables, evaluated at two different points. The two equations (1.1) 
define a difference equation, as well as a lattice. The general solution of the scheme 
depends on two constants and has the form 

y = y{n,Ci,C2), x = Ci, C2). (1.3) 
Equation (1.3) can also be rewritten as: 

y = y{x,Ci,C2), X ^ x{n,Ci,C2). (1.4) 
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The formula for y can be interpreted as interpolating between the points of the lattice 
and determining the solution y for all values of x. 

A standard equally spaced lattice is given by a specific choice of one of the two 
equations (1.1), namely E2 = x+ — x — h = 0, where ft, is a constant (the lattice 
spacing) and its solution is x„ = nh + Xo (where xq is one of the two integration 
constants of the scheme). 

Instead of the variables x, y, x+,y+, we can use 

u y+-y .1 <x 

X, y, h = x+- X, = . (1.5) 

x+ — X 

The continuous limit ft ^ of the scheme (1.1) is then obvious, namely, one of the 
equations should turn into a first order ordinary differential equation (ODE), the other 
into an identity (Uke = 0). We shall write the ODE as 

Eix,y,y')^y' -F{x,y)^0. (1.6) 

The symmetries and solutions of the difference scheme (1.1) that we shall obtain 
below will in the Unfit ft ^ turn into Lie point symmetries and solutions of the ODE 
(1 .6). To establish the connection, let us recall some well known results on symmetries 
of first order ODEs [1]. 

The Lie point symmetry group of the ODE (1.6) is always infinite-dimensional [1]. 
Its Lie algebra, the "syimnetry algebra", is reaUzed by vector fields of the form 

X = ^{x,y)d^ + <l>{x,y)dy (1.7) 

satisfying 

pvX{E)\^^^ = 0. (1.8) 
In (1.8) pvX is the first prolongation of X, i.e., [1] 

prX = ax,y)d:, + Hx,y)dy + r{x,y,y')dy,, (1.9) 
= cl>, + {cl>y-^,)y'-^y{y'f, 

where the subscripts are partial derivatives. Equation (1.8) amounts to a single first 
order linear partial differential equation for the two functions and (f) and as such it has 
infinitely many solutions. This is the reason why the Lie point symmetry of equation 
(1.6) is infinite-dimensional [1]. In general it may be difficult or impossible to find any 
explicit analytical solution; as difficult as finding an integrating multiplier. However, if 
we find at least one particular explicit solution of equation (1.8), we can obtain a one- 
dimensional subalgebra of the synnmetry algebra of equation (1.6). This is sufficient to 
integrate equation (1.6) in quadratures. 

All elementary methods of solving first order ODEs amount to special cases of the 
above procedure. 

A different application of the vector field X and its prolongation (1 .9) is to construct 
first order ODEs that are invariant under a given Lie group of local point transforma- 
tions, namely, those generated by the vector field X. In this case the functions ^{x, y). 
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(j){x. y) and hence also (f>^ (x, y. y'), are known. The invariant equation is obtained by 
solving the first order partial differential equation: 

Ux, y)d, + cj>{x, y)dy + cl)%x, y, y')dy,]E{x, y, y') = (1.10) 

for the function E{x, y, y'). Solving by the method of characteristics, we obtain two 
elementary invariants: 

dl 

h = h{x,y), h = h{x,y,y'), ^7^0. (1.11) 

dy' 

An invariant equation is given by any relation between Ji and I2, i.e., 

E{h,l2)=0, TTT^O- (1-12) 

If we are given a two-dimensional Lie algebra of vector fields {Xi,X2} and we 
require invariance under the two-dimensional Lie group they generate, then two differ- 
ent possibilities can occur. The first is that the two equations (1.10) (for Xi and X2 
respectively) have a common solution 

dl 

I = I{x,y,y'), g^T^O- (1-13) 

The invariant ODE then is 

F{I) = 0, (1.14) 

where F is arbitrary. We say that equation (1.14) is "strongly invariant" with respect to 
the group generated by {Xi,X2}- If no such invariant I{x, y, y') exists, then we look 
for an invariant manifold and a "weakly invariant" equation. This is obtained from the 
condition that the two equations ptXi{E) = 0, pvX2{E) = should be equivalent. 
This is a condition on the matrix of coefficients, i.e., 

rankf f ) = 1- d-lS) 

This condition, together with, say, prXi (E) = 0, provides the weakly invariant equa- 
tion. 

In Section 2 we will adapt the above results to the case of the difference system 
(1.1) and in Sections 3-9 consider many examples. The examples will be difference 
analogs of ODEs with known symmetry groups (linear equations, separable equations, 
etc). 

In each case we first present a one- or two-dimensional symmetry algebra of the 
ODE and use it to solve the equation. We have not found these symmetry algebras in 
the literature, but they are implicit in standard integration procedures (we have made 
them explicit). In each case we spell out the invariant difference schemes and choose 
one that has exactly the same general solution as the ODE. 
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2 Lie point symmetries and first order difference schemes 

The point of view that we will be taking here is the same as in previous pubUcations, 
e.g. [2-12] and references therein. Namely, Lie point symmetries of difference equa- 
tions will be continuous point transformations {x, y) (a;, y), taking solutions of the 
system (1.1) into solutions. They will be induced by a Lie algebra of vector fields of 
the same form (1.7) as for differential equations. The prolongation of the vector field 
will be different. No derivatives figure in equation (1.1); instead we prolong to other 
points on the lattice. In the case of the system (1.1) we have 

X = ^{x,y)d:jo + (j){x,y)dy + £,{x+,y+)d^^ + (l){x+,y+)dy^. (2.1) 

The continuous hmit (1 .9) is recovered by putting 

x+=x + h, y+ =y{x+) = y{x) + hy'{x) -\ (2.2) 

then, expanding into Taylor series 

£.{x+,y+) = £,{x,y) + h^a:{x,y) + hy'^y{x,y) ^ , (2.3) 

(t>{x+,y+) = <t>(x, y) + h(j)^{x, y) + hy'(t>y(x, y) -\ 

and acting with pr^X on a function of x, y, h and y^ (see equation (1.5)). We obtain 

pi^X{F{x, y, h, y^)) = {^(.t, y)dx + (pix, y)dy + [(j),^ + {(j)y - ^^)yx 

-Uyxf]dyjF{x,y,h,y,)+0{h) (2.4) 

Thus we have 

lim pi^X = vtX. (2.5) 

as required. 

Using the prescription 

pr^X(K)|^^^^^^^ = (2.6) 

we can find the symmetries of a given system (1.1). Instead, we shall start from an 
ODE and its known symmetries and construct the invariant difference scheme from the 
known vector fields. 

For a one-dimensional symmetry algebra we find invariants by solving the partial 
differential equation : 

[^{x,y)dx + (p{x,y)dy + ^{x+,y+)da;+ + (l>{x+,y+)dy_^_]F{x,y,x+,y+) = (2.7) 

by the method of characteristics. The elementary invariants are: 

I?=I?{x,y), =I^{x+,y+), =I^{x,y,h,y,), (2.8) 

with h and yx as in equation (1.5). Any two relations between the expressions (2.8) 
will give an invariant difference scheme, for instance 

I^=F{IP), 1^ = 1^ (2.9) 
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with F chosen to obtain the correct continuous Umit (we will drop the superscript D 
below). 

If we have a two-dimensional symmetry algebra, we will obtain two invariants. If 

one of them is of the type in (2.8), we again obtain a (strongly) invariant difference 
scheme. If not, we must look for an invariant manifold, given in this case by the rank 
condition 

rankf ^i^^i M^+.y+j Mx+,y+) )^^_ (2.IO) 

V Mx,y) Mx+,y+) J 

We will see below for specific examples that for each invariant ODE we obtain 
invariant difference schemes with the same invariance group and the same general so- 
lution. Indeed, our choice of the invariant difference schemes 

Ea{Ii,l2,l3)=0, a =1,2 (2.11) 
will be guided by two considerations: 

i) To obtain the original ODE in the continuous limit 

ii) To obtain a difference scheme that has exactly the same general solution as the 

original ODE (for any value of the lattice spacing h = x+ — x). 



3 Linear equations 

Let us consider the first order linear inhomogeneous ODE 

y' = a(x)y + b{x). (3.1) 

For convenience we redefine the given functions a{x) and b{x), putting a{x) = A'{x), 
b{x) = B'{x)e"^^''\ Equation (3.1) and its general solution are then written as 

y' = A'{x)y + B'{x)e^^-\ (3.2) 
y{x) = (B(a;) + fc)e^(-), (3.3) 

where k is the integration constant and the primes indicate a;-derivatives. 

Equation (3.2) has a two-dimensional Lie point synnmetry group, generated by the 
vector fields 

Xi = e^^^'^dy, X2=[y- B(x)e^(")] dy. (3.4) 

Now let us look for a difference scheme invariant under the group generated by the 
Lie algebra (3.4). The prolongations of Xi and X2 to the space {x, y, x+,y+} have 
only two invariants x and x^. However they do allow an invariant manifold, given by 
the condition (2.10), which in this case reduces to 

y^e-M^+) _ y Q-A(x) _ ^^^^^ ^ ^ Q (3 5) 

Adding an invariant lattice equation, e.g., 

x+ — X = h (3.6) 
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we obtain an invariant difference scheme (3.5), (3.6). Not only does this system reduce 
to the ODE (3.2) in the continuous limit, but equation (3.3) gives the exact general 
solution of the discrete system (3.5), (3.6) (for any value of h). 

Notice that while equation (3.5) is linear in y, it is not the difference equation one 
would get by the usual numerical discretization. To compare the two, let us reintroduce 
the discrete variable n, putting x = Xn, a;+ = Xn+i, y = yn> y+ = 

A "naive" discretization would be: 



y-+'_y- = A'ixr^)yn + S'(a.„)e^(-) (3.7) 



or possibly 



yn+l - Vn _ A{Xn+l) -A{Xn) ^^ B{Xn+l) - B{Xn) ^A(x„) 

— i/n T e l,J-oJ 

Xii-\-l Xji-^i Xfi 

with Xn = nh in both cases. The system (3.5), (3.6), on the other hand, is rewritten as; 

y„+ie-^(^"+i) - yn e-^^^-^ - B{xn+i) + B{xn) = (3.9) 

All of the above discretizations coincide in the limit h 0, however only the 
discretization (3.5), (3.6), i.e. (3.9), is exact in the sense that it has exactly the same 
general solution (3.3) as the ODE (3.1). 

Similar comments hold for all the discretizations that we present below in Sections 
4-9. We shall not repeat them each time. 

4 Separable equations 

Let us consider the separable ODE 

y' = f{x)g{y) (4.1) 

and for convenience redefine f{x) = A'{x), g{y) = 1/B{y) (the prime is an x- 
derivative, the dot a y-derivative). The (implicit) general solution of equation (4.1) 
is 

B{y) = A{x) + k, (4.2) 
where fc is a constant. The equation itself is rewritten as: 

y' = ^. (4.3) 
B{y) 

Equation (4.3) has a two-dimensional symmetry group, generated by 

Xr = ^.., X. = J^f.. i^A) 
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The group invariants of the group generated by the algebra (4.4) in the discrete 

space are 

h = A{x+) - A{x), h = B{y+)-B{y). (4.5) 
Using 7i and I2, we write an invariant difference scheme as 

B{y+) - B{y) - A{x+) + A{x) = 0, (4.6) 

A{x+)-A{x) = e{h), (4.7) 

where e is some constant, satisfying e(/i) for h ^ 0. Equation (4.2) clearly 
provides the general solution of equation (4.6). 
As a specific example, let us choose 

B{y) = y^, A{x)=x''. 

The difference scheme 

x^-x^ = e (4.9) 

is solved by 

y=(a;^ + fc)V^, a;„ = (ne + (4.10) 

where we choose a > 0, xo > 0. Clearly y{x) as in (4.10) also solves the ODE 
obtained in the continuous limit, namely 



y=T7-M^- (4.11) 



5 Exact equations 

We wiU consider in this section exact equations, that is, equations of the form: 

= i.e., A{x,y)dx + B{x,y)dy = (5.1) 

satisfying 

Ay = B^, i.e., A{x,y) =Va:{x,y), B{x,y) = Vy{x,y) (5.2) 

for some function V{x, y). 

Equation (5.1) is invariant under a one-dimensional group generated by 

X = B{x,y)d,-A{x,y)dy. (5.3) 

The general solution of (5.1) is given implicitly by the relation 

V{x,y) = k, (5.4) 



7 



where k is an integration constant. 

In the discrete case equation (2.7) leads to the characteristic system 

da; dy dx+ dy+ 



(5.5) 



Vy{x,y) Voo{x,y) Vy^{x+,y+) V^^{x+,y+) 
and hence to three invariants 

h=V{x,y), h = V{x+,y+) (5.6) 

^ J Vy_^{x+,y+{x+,l2)) J Vy{x,y{x,Ii))' 

To obtain the integrals involved in I3 we have solved the equations (5.6) for y and y+. 
The discrete version of equation (5.1) is: 

V{x+,y+)-V{x,y)=0 (5.8) 

with (5.4) as its solution. Equation (5.7) can be used to define the invariant lattice. 
As a specific example, consider the ODE 

y' = ^, \ . - 1- (5.9) 

2{x + y) 

It is exact and we have 

V ={x + yf -X. (5.10) 
The vector field (5.3) in this case is 

X = 2{x + y)d^-{2{x + y)-l)dy. (5.11) 

The invariants in the discrete case are 

h = {x + yf -X, h = {x+ + y+f - x+, (5.12) 
h = x+ + y+ - X - y. 

The invariant difference scheme can be written as 

{x++y+f -(x + yf ^x+-x, (5.13) 
x++y+- x-y = e. (5.14) 

Returning to the general case, we see that the discrete analogue of an exact ODE is 
equation (5.8). The invariant lattice can be given by equation (5.7). In the continuous 
limit, we have /s — > 0, i.e., Xj^ ^ x and (5.8) goes to 

V,+Vyy, = Q. (5.15) 
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6 Homogeneous equations 

Let us consider the first order ODE 



2/' = x'=-1f(^), (6.1) 

wliere F is an arbitrary smootli function and A: is a real constant. This is the most 
general first order ODE invariant under the scaling group 

y = e^'^y, X = e^x (6.2) 

generated by the vector field 

X = xdx + kydy. (6.3) 

For convenience, we replace the function F{t) in equation (6.1) by F{t) = 1/H{t) + 
kt, so equation (6.1) is rewritten as 

y> = ^!ll^+kl. (6.4) 

The change of variables {x, y{x)) {t, z{t)) with 

t=^, z = logx, a; = e^ y = te'^^ (6.5) 
x" 

will straighten out the vector field X and transform equation (6.1) into 

zt=H{t). (6.6) 

Solving equation (6.6) and returning to the original variables, we obtain the general 
solution of equation (6.1) in the form 

y{x) = x^H-\\ogx-C), (6.7) 

where C is an integration constant and H^^ is the function inverse to H{t). 

Let us now find the invariant difference scheme corresponding to equation (6.1). 
Prolonging the vector field X of equation (6.3) as in equation (2.1), we find three 
elementary invariants 

h = —r, h = —r, h = — • (6.8) 

Using them we write an invariant difference scheme as: 

\ogh-H{h) + H{h) = 0, (6.9) 
/g-l-e = 0. (6.10) 

More explicitly, we have 

logx+-logx-ij(^)+i?(^) = 0, (6.11) 

Xj^—x — ex = 0, (6.12) 
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where e is some constant. The e ^ Hmit of the difference scheme (6.11), (6.12) is 
the ODE (6.4), as required, and its general solution is (6.7), together with 

Xn = i€ + irXo. (6.13) 

As in the previous examples, the invariant difference scheme has the same exact solu- 
tion as the original ODE. 

As a specific example, we take k = 1 and the ODE 

y' = (6.14) 

xy 

The difference scheme in this case is 

M'-('-y=2log^-±, ^ = e. (6.15) 
\x+ J \x) X X 

The solution of both (6.14) and (6.15) is 

y = x^1\ogx-G (6.16) 
(together with (6.13) in the discrete case). 

7 Rotationally invariant equations 

Let us consider a rotation in the (x, y) space. It is generated by 

X = ydx-xdy. {lA) 
The most general first order ODE invariant under these rotations is: 

y'-'^TWT^ P=^^^. (7.2) 

y + K (p)x 

To solve equation (7.2) we straighten out the vector field (7.1) by going to polar coor- 
dinates {x,y{x)) {p,a{p)) 

x = pcosa, y = psina. (7.3) 

Equation (7.2) reduces to 



and we obtain 



For simplicity, let us restrict to the special case K{p) = K = constant. We then 
obtain the solution of (7.2) as a logarithmic spiral 

Ka -h log p = log po (7.6) 
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or, in the original variables 

1 y 

- \og{x^ +y^)+ K arctan - = log po- (7.7) 
Now let us consider the discrete case. The prolonged vector field is: 

prX = ydx - xdy + y+d^^ - x+dy^ . (7.8) 

This means that the pairs (a;, y) and {x+, y+) will transform like vectors undergoing a 
rotation in a Euclidean plane. We can form 4 invariants 

Ii=x'^+y'^, I2 = x1+yl, h=xy+-x+y, h = xx++yy+ (7.9) 

with one relation between them, namely 

Il + Il = hh. (7.10) 

As an invariant difference scheme, we write 

El = log + 2i4: arctan = 0, (7.11) 
h h 

E2 = ^{l2-Ii)+Kh = (7.12) 

with K a constant. Using equation (2.2) we can check that the continuous limit of both 
these expressions is the ODE (7.2) (with K constant). In general we have 



d{x+,y^ 



^ 0, lim 



d{Ei,E2) 



d{x+,y^ 



0. (7.13) 



Expression (7.7), the exact solution of the ODE (7.2), is also an exact solution of 
the system (7.11), (7.12). 

The situation is more transparent in polar coordinates (7.3). For K = constant, the 
ODE is d 1 

dp Kp ^ ^ 

The prolongation of the vector field corresponding to rotational invariance in the dis- 
crete case is 

prX = a„ + S„+. (7.15) 

The invariants are 

h=a+-a, l2 = p+, h=P- (7.16) 
An invariant difference scheme corresponding to the ODE (7. 14) is 

a+-a = -^{\ogp+-\ogp), (7.17) 

p+-p = e. (7.18) 

The exact solution of (7.14) and (7.17) is given by equation (7.6). The lattice deter- 
mined by equation (7.18) is uniform: 

p„=en + pQ. (7.19) 
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8 Invariant difference schemes on uniform lattices 

Let us consider the ODE 

y' = -^ + f'i^)4-^ ^v^Q. (8-1) 

where A{x, y) and f{x) are some smooth functions. Equation (8.1) is invariant under 
transformations generated by 



1 

Ay{x,yy 



X=— -dy. (8.2) 



As a matter of fact, any first order ODE invariant under transformations generated by 
a vector field of the form 

X = (j>{x,y)dy (8.3) 

can be reduced to the form (8.1). 

The general solution of equation (8.1) can be written imphcitly as 

A{x,y) = f{x) + C. (8.4) 

In the discrete case the first prolongation of the vector field (8.2) has 3 invariants; 
they can be chosen to be: 

h=x+, l2 = x, l3 = A{x+,y+)- A{x,y). (8.5) 

An invariant difference scheme on a uniform lattice, having the ODE as a continuous 
limit is: 

A{x+,y+) - A{x, y) = f{x+) - f{x), (8.6) 

x+ — X = h. (8.7) 

Its solution is given by equation (8.4) together with 

Xn = hn + xo- (8.8) 
As a specific example, let us choose 

X = xydy. (8.9) 

The invariant ODE in this case is 

y' =y-\ogy + q'{x)xy (8.10) 

X 

the corresponding difference scheme is 

— ^ogy+--\ogy-q{x+)+q{x)=Q, x+ - x = h. (8.11) 

x+ x 



12 



The solution of both (8.10) and (8.11) is 

y = ex(c+9(x))^ (8.12) 

where C is an integration constant. 

The vector field (8.2) is not the only one compatible with a uniform lattice. Another 
one is 

X = + (t){x,y)dy (8.13) 
for any function (j){x, y). As an example, let us take 

X = 9,+a;V5j/- (8.14) 
The corresponding invariant first order ODE is 

y' = {k{0+x'^)y\ Q=-L_x'^+^ + ^y-'+\ (8.15) 

where k is any function of For k = ko =constant. The solution of equation (8.15) is 

_ 1 

1 



y={l-b)-- 



kox + 



a+1 



(8.16) 



where C is an integration constant. 

The invariants of the discrete prolongation of the vector field (8.14) are 

^a+l -6+1 -6+1 

h = x+-x, I^ = —- + !L-—^ /3^^ + ^i_. (8.17) 

a+1 0—1 a + 1 0—1 

We write an invariant difference scheme as 

h-h- koh =0, h = h (8.18) 
with ko and h constant. More explicitly (8.18) is 

Yz-b^yl-' - y'-'] - Y^K' - ^"^'] - ko{x+ ~x)^o, (8.i9) 

x+-x = h. (8.20) 

The continuous limit of (8.19) is the ODE (8.15). Moreover, equation (8.16) is the 
exact solution of (8.19) for any value of h (not just the limit h —>■ 0). 

9 Invariant difference schemes on exponential lattices 

The vector field (8.2) is also compatible with an invariant scheme of the form (8.6) on 
the exponential lattice 

^^^=e, i.e. Xn = {e + irxo (9.1) 
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with e = constant. Taking the Hmit e ^ we again obtain the ODE (8.1). 

Another synunetry, leaving the lattice (9.1) invariant, is generated by the vector 
field 

X = xd^+<i){x,y)dy (9.2) 

with (/)(a;, y) arbitrary. 

As an example, consider 



X = xd^ + x^y^dy, a ^ 0, 6 ^ 1. 
The corresponding invariant ODE is 

1 



1 



1 



X a — 1 



For k = ko constant, the solution is 



y={l-b)- 



/co log cc H \-C 

a 



In the discrete case, the prolongation of X has three invariants 



T T 

-'1 ~ ' -'2 I, 1 

X 0—1 



yl-h ,j^a 



h = 



b-1 



J_ 
a 



(9.3) 



(9.4) 



(9.5) 



(9.6) 



An invariant scheme having (9.4) with fc = fco as its limit and (9.5) as its solution is 

h-h- ko log/i =0, h-l^e, (9.7) 

i.e.. 



y 



l-b 



y^ 

1-6 1-6 



X 

a 



A;o(loga;+ - logo;) = 0, 
x+ — X = ex+. 



(9.8) 
(9.9) 



10 Conclusions 

Looking directly for synmietries of a difference scheme of the type (1.1) is not a par- 
ticular fruitful enterprise. Just as in the case of first order ODEs one gets an underde- 
termined system of equations. Infinitely many solutions exist, but there is no algorithm 
for finding them. 

We have taken the complementary point of view. We have postulated the form 
of a vector field, then found ODEs and difference schemes, invariant under the corre- 
sponding symmetry group. The symmetry makes it possible to solve the ODE exactly 
analytically. The invariant differential scheme then has the same solution. More pre- 
cisely, the symmetry leads to a family of difference schemes, one of which has solutions 
coinciding with those of its continuous limit (the original ODE). 
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Essentially we have constructed a partial "catalogue" of exactly solvable two-point 
schemes. This corresponds to a list of exactly solvable first order ODEs: linear equa- 
tions, separable equations, exact equations, homogeneous equations, etc. The "com- 
plete list" is infinite: for any chosen realization of a one- or two-dimensional Lie alge- 
bra we can construct an invariant ODE and an invariant difference scheme. 

Without the symmetries to guide us, the obtained difference scheme are not obvious 
at all. Starting from an ODE (1.6) and discretizing in a standard way, i.e., replacing the 
derivative dy/dx by a discrete derivative and writing 

'^"^'^ = F{xn,yn), Xn = hn + xo (10.1) 

we would loose virtually all symmetries. Moreover, the exact solution of equation 
(10.1) would differ from that of the ODE (1.6) by terms of the order h. For our sym- 
metry dictated difference schemes, the exact solutions coincide with those of the ODEs. 
The main result of this article can be summed up as follows. 

Theorem. For every first order ODE there exists an invariant two point difference 
scheme with exactly the same general solution as the ODE. 

Proof. Consider the ODE (1.6) and assume that we know its general solution in the 
form of a first integral 

h{x,y) = K, hyy^O. (10.2) 

The general element (1.7) of the symmetry algebra of equation (1.6) must annihilate 
the function h{x, y). The equation Xh = implies that X has the form: 

X = ax,y){d,~^^d,)i, (10.3) 

where ^(x, y) is an arbitrary smooth function. 

Let us now find the invariants of the group action induced by the vector field (10.3) 
in the space (a;, a;+, y, y+). They are obtained by solving the characteristic system: 

da; _ da;+ _ hy{x,y)dy _ hy_^_{x+,y+)dy+ ^^^^^ 



^{x,y) ^{x+,y+) -£,{x,y)K{x,y) -^{x+,y+)K^{x+,y+)' 
The invariants hence are 

h = h{x,y), l2 = h{x+,y+), 

J £,(x,y(x,h)) J £,{x+,y+{x+,l2)) 

where the function ^(.t, y) can be freely chosen. 

A difference scheme that has (10.2) as its general solution is 

/i = I2, I3 = c, (c = constant). (10.6) 

QED 

Three comments are in order here. 
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1. The proof given above is not a constructive one. It assumes that we already 
know the general solution (10.2) of the ODE. This was not assumed in the rest 
of this article, where we constructed the difference schemes using only one, or 
sometimes two elements of the symmetry algebra of an ODE. 

2. The arbitrariness in the function ^(.t, y) can be put to good use in the choice of 
lattices. For instance, choosing ^{x, y) = 1 we obtain a uniform lattice as in 
(8.8). Choosing ^ = x we obtain an exponential lattice, as in equation (9.1). 

3. The theorem and the results of this article are specific of first order ODEs and 
two point difference schemes. 

Let us compare with the case of second order ODEs and three-point difference 
schemes. S. Lie classified all (complex) second order ODEs into equivalence classes 
according to their Lie point symmetries [13]. A similar classification of three-point 
difference schemes is much more recent [11]. It was shown [12] that if a three-point 
difference scheme has a symmetry group of dimension 3 (or larger) with at least a two- 
dimensional subalgebra of Lagrangian symmetries, then the scheme can be integrated 
analytically. A crucial element in the integration was the existence of a Lagrangian and 
the interpretation of the difference scheme as a discrete analogue of an Euler-Lagrange 
equation [11, 12]. 

The fact that we obtain differential equations and difference schemes that have 
identical symmetries and solutions has interesting implications. It suggests a certain 
duality between continuous and discrete descriptions of physical phenomena. Thus, 
exactly the same physical predictions may be described by some continuous curve, or 
by a series of points on this curve, distributed with an arbitrary density. 
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